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Abstract 

We report temperature-dependent surface x-ray scattering studies of the orbital ordered surface 
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in LaosSri 5MHO4. We find that the interfacial width of the electronic order grows as the bulk 

ordering temperature is approached from below, though the bulk correlation length remains un- 

changed. Close to the transition, the surface is so rough that there is no well-defined electronic 

j> surface, despite the presence of bulk electronic order, that is the electronic surface has melted. 

00 

Above the bulk transition, finite-sized isotropic fluctuations of orbital order are observed, with 
a correlation length equal to that of the electronic surfaces' in-plane correlation length at the 
Q transition temperature. 

PACS numbers: 73.20.-r 
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Bulk crystals can disorder on warming through one of two scenarios: surface melting or 
surface freezing. In the first, the surface of the solid melts at a lower temperature than the 
bulk. This is a common and even useful property of condensed matter and is, for example, 
what allows a skier to glide across the snowjTJ.In the second, the bulk melts at a lower 
temperature than the surface [21 EJj. An as yet unanswered question is what happens when 
the order is not that of a crystal lattice, but is rather that of electronic order? We show, in 
this letter, that orbital order in the manganite Lao.5Sr 15 Mn04 is an example of electronic 
surface melting. 

The temperature dependence of electronic order at a surface or interface is of fundamen- 
tal interest since it can provide important insights into the nature of electronic interactions, 
and in particular the relative strength of the surface and bulk interactions. Nor is the ques- 
tion entirely academic. Next-generation devices that rely on electronic order will require 
understanding and control of electron transport across such surfacess at finite tempera- 
ture, and this interfacial behavior will also have large implications for the properties of any 
nanoscaled strongly correlated materials at finite temperature. Thus this is a question of 
some importance. 

Here, we use orbital order at the surface of a strongly correlated transition metal oxide 
as an exemplar system for studying electronic order at a surface. Before proceding, it is 
important, to define what we mean by a surface. The crystallographic surface is the surface 
between the solid crystal and the vacuum, as obtained by cleaving the crystal. The electronic 
surface (in this case the orbital surface) is the surface between the orbitally-ordered bulk 
and an orbitally disordered region located below the crystallographic surface. There has 
been much work associated with the ground state behavior of such surfaces, starting with 
the original suggestion that "electronic reconstruction" can occur at the electronic order- 
disorder interface [4] . Surface effects in electronic order have been seen in such cases as 
the ferromagnetic surface of bilayer manganitesjS], the interface between ferromagnetic and 
superconducting oxides [6], the orbital surface of manganites[7] and the interface between 
LAO/STO heterostructuresjE]. However, to date there has been little work on the evolution 
of the electronic surface on warming through the bulk ordering temperature. The purpose 
of this work is to redress this imbalance. 

We report temperature-dependent surface x-ray scattering studies of the orbitally ordered 
surface in Lao.5Sr 15 Mn04 (LSMO). We find that the interfacial width of the electronic order 
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grows as the bulk ordering temperature is approached from below, while the bulk correlation 
length remains unchanged, suggesting that the surface starts to roughen while the bulk is 
well ordered. Close to, but below, the transition, the surface is so rough that there is no 
well-defined electronic surface despite the presence of bulk electronic order: the electronic 
surface has melted. Above the bulk transition, finite-sized isotropic fluctuations of orbital 
order are observed, with a correlation length equal to that of the surface in-plane correlation 
length at the transition temperature. These results have important implications for the 
behavior of electronic surface at finite temperatures and we shall discuss these in view of 
potential future devices. 

LSMO has been studied extensively in the bulk, revealing long range orbital order below 
~240 K, with 3D magnetic order appearing below ~120 K (9J [TO]. This orbital ordering 
gives rise to concomitant lattice distortions and thus may be studied utilizing x-ray scat- 
tering techniques. Specifically, by measuring the distribution of the scattering around the 
orbital-order superlattice reflection, one obtains information about the bulk order param- 
eter. Information about the orbital surface comes from the "orbital truncation rods" [7] . 
These are rods of scattering normal to the surface connecting orbital superlattice peaks and 
are analogous to crystal truncation rods [HI [12] which result from the termination of a 
3D crystal by a surface. The intensity distribution of the orbital-truncation-rods provides 
quantitative information on the interfacial roughness (interfacial width) and the in-plane 
correlations at the orbital surface. 

The experiments were carried out on beamline 6ID-B at the Advanced Photon Source. 
The incident photon energy was 12 keV. A single crystal of Lao.sSri^MnC^, grown in a 
floating zone furnace, was post-cleaved in air to reveal a [001] surface and immediately 
placed into the vacuum of a closed-cycle refrigerator. Q z is along the surface normal. Scans 
were performed using a 4+2 diffractometer[13], accessing to the (2.25, -0.25, L) orbital rod 
and the (2.25, -0.25, 2) orbital Bragg peak. An incidence angle of 0.6 ° was used for all 
measurements. The sample was indexed with a = b = 3.78 A and c = 12.4 A. Momentum 
transfer is either denoted as Q x , Q z in units of A -1 , or L, in units of r.l.u., where 1 r.l.u. 
= 0.507 A" 1 . Measurements of the (2.25, -0.25, L) orbital-truncation rod were made by 
scanning along Q x , perpendicular to the Q 2 direction. A 2-D mesh of the intensity obtained 
from a series of such scans, collected at 170 K, is shown in the top right panel of Fig. [TJ The 
top left hand side of Fig. [T] shows a simulation of the expected intensity of the superlattice 
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FIG. 1: (color online) (Top left) Simulated x-ray scattering in the (Q x ,Qz) plane for a 2D 
Lorentzian-squared function, with widths equal to those of the bulk orbital order reflection. (Bottom 
left) Line cuts at Q 2 = A -1 (blue) and Q 2 = 0.28 A -1 (red) through the simulated x-ray scat- 
tering intensity. (Top right) 2D Intensity plot of the measured intensity for the orbital truncation 
rod in Lao.5Sri.5Mn04 at 170 K. The streak of scattering, elongated along L due to the truncation 
of the orbital order by its surface can clearly be seen. (Bottom right) Line cut at Q z = 0.28 A -1 
through the orbital truncation rod. 
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reflection, in the absence of the orbital truncation rod, constructed from a 2D Lorentzian- 
squared function with Q z and Q x widths equal to the values measured at the (2.25, -0.25, 
2) bulk superlattice reflection. The bottom panels show line cuts through the respective 
meshes. It is immediately apparent that the two behaviors are quite different; the observed 
scattering does not originate from tails of the orbital Bragg peak, but rather is the orbital 
truncation rod[7j. 

Figure [2] shows two Q x . scans through this orbital truncation rod, taken at Q z = 0.1 A -1 
and 0.28 A -1 , for two temperatures, T=170.0 K and 223.5 K, i.e. well below, and just 
below, the bulk ordering temperature, Tqo = 226 K. At low temperatures, a two-component 
lineshape is observed in the truncation rod scattering. This immediately implies that the 
orbital surface has a well-defined average height but finite in-plane correlations at the sur- 
face [TT]. We refer to these two components as the sharp and broad components. In 
the theory of a nearly-smooth surface [11], the sharp component results from the average 
surface, and the intensity variation of this component as a function of Q z is a measure of 
the roughness, or interfacial width, of that surface. A rough surface results in a faster decay 
of the sharp component intensity with Q 2 , compared to that of a perfectly smooth surface. 
The broad component on the other hand, results from the finite in-plane correlation length 
of the surface, with a width inversely proportional to that correlation length. The rate of 
increase of the integrated intensity of the broad component as a function of Q z increases 
with increasing interfacial roughness. We find that the two length-scale scattering is best 
fit by a Lorentzian-squared function for the sharp component and a Lorentizan function for 
the broad component: 
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where Is and 1b are the peak intensities of the sharp and broad components respectively, 
and Ts and Tb are the respective widths. Correlation lengths £, are defined as £ = 1/r. 
For a perfect crystal, the presence of an average surface would give rise to a delta-function 
sharp component [TT]. For the present case, with a finite bulk correlation length, we model 
the sharp component with a Lorentzian-squared function [7] and constrain the width of this 
component to be equal to the width of the bulk orbital order superlattice reflection, Vg= 
0.0055 A -1 . 

Taking a close look at the T=170 K data, we see that at large Q 2 , the broad component 
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FIG. 2: (color online) (Top panel) Q x scans through the orbital truncation rod at 170 K and 223.5 

K (red circles and blue squares respectively). Data taken at Q z = 0.1 A -1 (Top panel) and = 

0.28 A -1 (Bottom panel). The results of the fits to Eq. 1 are shown as solid lines. 
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has a much greater integrated intensity than the sharp component, while at small Q^, the 
sharp component dominates. This is the expected behavior for a nearly smooth surface. 
Note that while the broad component is much broader than the sharp component, it is not 
as broad as would be expected if this was a line cut through the tails of the orbital Bragg 
peak (bottom panel, Fig. [I]), that is both the sharp and broad components are associated 
with the orbital truncation rod, and hence reflect the properties of the orbital surface. 
At low temperatures, we find that the in-plane correlation length of the orbital surface is 
1/Tb=13 A, significantly shorter than the bulk correlation length (1/T S =180 A), and in 
agreement with earlier work [7]. 

Having established that we are probing the orbital surface, we can now study the temper- 
ature dependence of that surface upon warming through the bulk orbital order transition. 
The T = 223.5 K data of Fig. [2] already reveal dramatic changes from the low temperature 
lineshape. Specifically, there is no longer any evidence of a sharp component at any Q z , and 
the broad component is both broader and more intense than it was at low temperatures. 

Figure [3] shows the evolution of this behavior in detail. In the top panel, the temperature 
dependence of the integrated intensity of the bulk superlattice Bragg peak (L = 2), the 
sharp component of the surface scattering, measured at = 0.13 A -1 , and the broad 
component as measured at Q z = 0.2 A -1 , are shown. As the bulk orbital ordering transition 
is approached, the amplitude of the bulk superlattice peak decreases, disappearing at T Q o = 
226 K. The integrated intensity of the sharp component falls faster than the bulk orbital order 
Bragg peak and becomes indistinguishable from zero at T=222 K, i.e below Too- I n contrast, 
the integrated intensity of the broad component smoothly increases. The correlation length 
of the bulk orbital order is also shown, measured along the and Q x directions, in the 
bottom panel of Fig. |3j along with the correlation length as measured on the the broad 
component, along the Q x direction. The latter is also plotted on a linear scale in the 
middle panel of Fig. [3] These data show that at low temperatures, the bulk orbital order is 
anisotropic, with the larger correlation length in the Q x direction. This correlation length is 
larger than that of the broad component, indicative of a larger in-plane correlation length 
in the bulk compared to the orbital surface. As the temperature is raised, the correlation 
length of the bulk orbital reflection remains constant until Too a t which point it rapidly 
decreases, becoming isotropic above the transition. In stark contrast, the corrleation length 
as measured on the broad component decreases steadily as the temperature is raised (middle 
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FIG. 3: (color online) (Top panel) Temperature dependence of the integrated intensity of the bulk 
orbital order superlattice reflection (red squares) measured at L = 2 , the integrated intensity of 
the sharp component of the orbital truncation rod measured at Q 2 = 0.13A -1 (blue circles) and the 
broad component measured at Q 2 = 0.2 A -1 (magenta diamonds). (Middle panel) Temperature 
dependence of the correlation length as measued on the the broad component of the orbital trun- 
cation rod at Q 2 = 0.2 A -1 . (Bottom panel) Temperature dependence of the correlation lengths 
of the bulk orbital order superlattice reflection measured along the Q 2 and Q x directions (green 
circles and red squares respectively) and the correlation length as measured at the on the broad 
component of the orbital truncation rod at Q z = 0.2 A -1 . 
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panel, Fig. [3J), indicative of a slow reduction of the in-plane interfacial correlation length. 
Above Tqo j the correlation lengths of the bulk orbital order Bragg reflection and the broad 
component become almost equal. 

Before discussing the significance of these experimental results, it is useful to first discuss 
what the possible behavior for the orbital surface might be, on approaching the bulk orbital 
order phase transition. One possibility is that of surface freezing, in which an ordered region 
exists at the surface of a disordered bulk. Soft matter systems such as liquid crystals [H] 
and alkanesj21 E] exhibit such behavior. A second, more common, scenario is that of sur- 
face melting (or pre-melting). Surface melting manifests itself as a disordered, liquid-like, 
region on the surface of a material, which grows in thickness upon approaching the bulk 
transition temperature from below [15]. In principle, electronic surfaces could exhibit either 
of these behaviors. In the study of magnetic surfaces, for example, there have been various 
experimental observations of differing surface and bulk magnetic phase transitions [161 E] • 
On the theoretical side, it was found that for the ferromagnetic Ising model, - a model that 
has been used for orbital order in the manganites|18j - it is possible to obtain either surface 
melting or surface freezing behavior, depending on the ratio of the exchange interactions for 
the surface and the bulk[19]. The question is then, which of these scenarios is the relevant 
one for the orbital surface discussed here? 

The data presented above reveal very different behavior for the electronic surface and 
the bulk orbital order on warming through the bulk orbital ordering transition: At T = 
170 K, well below the bulk transition temperature, the orbital surface is well defined but 
rougher than the crystal chemical surface, and has a shorter correlation length than the 
bulk orbital order. This situation is shown pictorially in Fig. |4]^a) where we represent the 
different correlation lengths by regions of blue and yellow shading. 

As the temperature is increased, the integrated intensity of the sharp component falls 
and the integrated intensity of the broad component increases, indicating that the interfacial 
width, or roughness, is increasing. Further, the correlation length measured on the broad 
component is found to steadily decrease, in contrast to the correlation length of the bulk 
orbital order, which remains constant throughout this temperature regime. This implies that 
the in-plane correlation length at the orbital surface is decreasing, while the bulk correlation 
length is unchanged. 

At a temperature just below the bulk orbital order phase transition, the integrated inten- 



9 



a) T « Too 

Crystallographic 
Surface 



Orbitally Ordered 
Surface 

Orbitally Ordered 
Bulk 

b) T = Too c) T » Too 

l l 



FIG. 4: (color online) Schematic representation of the crystalographic surface, the orbital surface 
and the orbitally ordered bulk for (a) T <C T OQ , (b) T « T OD and (c) T 3> T OD . The size of the blue 
and yellow regions in the orbitally orderd bulk represent schematically the correlation length of the 
orbital order in the bulk. The roughness of the orbital surface is shown as the height variations of 
the shaded (blue and yellow) surface and the in-place correlation length by the size of the shaded 
surface regions. The size of the all the colored regions are physically correct, both relative to each 
temperature snapshot and are consistent with the x-ray scattering data analysis. 
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sity of the sharp component has dropped below detectable levels, and the broad component 
completely dominates, indicating that there is no average orbital surface at this tempera- 
ture, and that the orbital surface has melted. This is shown pictorially in Fig. Rib). It is 
important to note that at this temperature, the bulk is still as well correlated as it was at 
low temperatures: Thus, this is a case of electronic surface melting. 

Increasing the temperature further, one crosses the bulk phase transition and the bulk 
correlation lengths decrease to have approximately the same magnitude as the interfacial 
correlation length. In the absence of an surface, this implies that the orbital order exists as 
fluctuations in a disordered background. Further increasing the temperature causes the size 
of these fluctuations to decrease along with the integrated intensity, until they are no longer 
observable. This case is shown in Fig. Qc). 

We conclude by noting that this observation of electronic surface melting has potentially 
profound implications on the physics of strongly correlated materials and their applications 
in any device. Our results show that, at least for this system, the electronic surface becomes 
increasingly rough and less well defined on warming, and vanishes below the bulk ordering 
temperature. Materials with electronic surfaces that behave in this way could be problematic 
in devices operating at finite temperature. For example, colossal magnetoresistance|20j is a 
potentially very useful property, offering a dramatic change in resistance in applied magnetic 
fields. However, in manganites, this effect requires the system to operate at a temperature 
close to the bulk transition temperature. The present work shows that this temperature 
corresponds to the roughest (thickest) orbital surface, which could be a problem for electron 
transport across that surface. Further, the properties of any nano-scaled device are likely to 
be dominated by the surface / interface behavior due to the large surface-to-volume ratio. 
Continuing with the CMR example, the stark difference between the surface and bulk at this 
temperature means that at the nanoscale, CMR may manifest in a very different form - if at 
all. Finally, we comment that while the data presented here applies only to the orbital surface 
of LSMO, the implications of this work are likely to extend much further. For example, we 
note that much electronic information about strongly correlated systems is obtained from 
inherently surface sensitive techniques such as STM and ARPES. For systems that exhibit 
behavior analogous to the present case of electronic surface melting, such techniques would 
measure the melted, rough layer, and not reveal the true bulk temperature evolution. Our 
results suggest that future experimental and theoretical investigations of the temperature 
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dependence of electronic surfaces are to be encouraged. 
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